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Abstract 

' We model the observed glueball mass spectrum in terms of energies for tightly knotted and 

I linked QCD flux tubes. The data is fit well with one parameter. We predict additional glueball 

■ masses. 

■ 

^ ' Introduction. — The interpretation of non-gg states is a puzzle with a long and controversial 

\ history pQ. Many experiments |2] report states that do not fit neatly into the quark model. These 

CO ■ states can be broadly classified as: (1) hybrids, which arc bound states of quarks and gluons like qqG 

\ with quantum numbers J^*^ — ^, 1 '", 1 , 2 . . .; (2) exotics, for example, four and six quark 

■ states such as qqqq and qqqqqq with quantum numbers J^'-^ = 0~~, 0"' , 1 ^, 2^ , . . .; (3) glueballs 
I with pointlike or collective (e.g., strings a la Niclscn-Olesen P], or flux tubes) glue. Glueballs do not 
' contain valence quarks, but there could be sea/virtual quarks within the glueball or in the currents 
I that support the flux tubes. From a bag model perspective one is led to suppose that the lightest 

^ • non-qq states are those with no constituent quarks, i.e., the glueballs. Lattice calculations, QCD 

I I sum rules, electric flux tube models, and constituent glue models leads to a consensus that the 

lightest non-gg states are glueballs with quantum numbers J++ = 0++ and 2++ @]. We will model 
(— I I all J++ states (i.e., all fj and fj states listed by the Particle Data Group (PDG) 0), some of which 

■ will be identified with rotational excitations, as knotted/linked chromoclcctric QCD flux tubes [S]. 
. !^ I Besides the fact they do not fit in the quark model [H] , glueballs have some other characteristic 

■ signatures, including: enhanced central production in gluon rich channels, branching fractions in- 
^ I compatible with qq decay, reduced 77 coupling, and OZf supprcsion. All the J^"*" states wc consider 

■ - - I have some or all of these properties. For instance, none have substantial branching fractions to 77. 

However, mixing with qq isoscalar states can obscure some of the properties. A number of candidates 
with masses below 2.5 GeV have been identified. Beyond their masses and widths, and some of their 
branching ratios much remains to be learned about these states. 

Knotted magnetic fields (which we will treat as solitons) have been suggested as candidates for 
a number of plasma phenomena in systems ranging from astrophysical, to atmospheric |Zj , to Bose- 
Einstein condensates [H| . The energies of these solitons are sometimes difficult to quantify since they 
depend on parameters of the plasma, including temperature, pressure, density, ionic content, etc.; 
however, we will argue that in QCD a well defined soliton energy can be identified. 

As has been shown in plasma physics, tight knots and links (defined below) correspond to 
metastablc minimum energy configurations. We will argue by analogy that quantized tightly knoted 
and linked QCD fiux tubes are glueballs. (In what follows, we often use the term "knots" to mean 
knots and/or links.) 

Movement of fiuids often exhibits topological properties (for a mathematical review see e.g. |H]). 
For conductive fiuids, interrelation between fiuid motion and magnetic fields via magnetohydrody- 
namics may cause magnetic fields, in their turn, to exhibit topological properties. For example, for a 
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perfectly conducting fluid, the (Abelian) magnetic helicity Lh = J d^x CijkAidjAk is an invariant of 
motion |10| , and this quantity can be interpreted in terms of knottedness of magnetic flux hnes . 

The dynamics of the magnetic flelds foUows the dynamics of the hquid (magnetic flux hnes 
arc "frozen" into the fluid), and one finds that a perfectly conducting, viscous and incompressible 
fluid relaxes to a state of magnetic equilibrium without a change in topology |12j . As a result, for 
topologically non-trivial plasma flows (with knotted streamlines), the "freezing" condition forces 
topological restrictions on possible changes in fleld configurations. For linked non-intersecting loops 
Ca with magnetic fluxes <i>Q, the hehcity becomes QT] Lh = X^a^^b L{Ca, Cbj'^a^b- Here L{Ca, Cb) 
is the Gauss linking number ^^j- By its topological nature, helicity in the QCD flux can be one of 
the quantum numbers characterizing giueballs. However, there is another invariant called the knot 
energy that is less obvious but as important in the classification of solitonic knots. 

Knot energies. — Consider a hadronic collision that produces some number of baryons and mesons 
plus a gluonic state in the form of a closed QCD flux tube (or a set of tubes) . From an initial state, 
the flelds in the flux tubes quickly relax to an equilibrium conflguration, which is topologically 
equivalent to the initial state. (We assume topological quantum numbers are conserved during this 
rapid process.) The relaxation proceeds through minimization of the field energy. Flux conservation 
and energy minimization force the flelds to be homogeneous across the tube cross sections. This 
process occurs via shrinking the tube length, and halts to form a "tight" knot or link. The radial 
scale will be set by ^qqjj • The energy of the final state depends only on the topology of the initial 

state and can be estimated as follows. An arbitrarily knotted tube of radius a and length / has the 
volume TTo^l. Using conservation of flux ^e, the energy becomes oc /(tr$|;) / (tto^). Fixing the radius 
of the tube (to be proportional to we find that the energy is proportional to the length /. 

The dimensionless ratio e{K) = I /{2a) is a topological invariant and the simplest definition of the 
"knot energy" [Tl] . 

Many knot energies have been calculated by Monte Carlo methods JS] and certain types can 
be calculated exactly (see below), while for other cases simple estimates can be made (see Table 
For example, the knot energy of the connected product of two knots Ki and K2 satisfies 

e(Ki#K2) <e(/^i)+£(i^2) (1) 

A rule of thumb is 

e{Ki#K2) ~ e{Ki) + eiK^) - (2^ - 4), (2) 

which results from removing two half tori, one from each knot, and replacing these with two con- 
necting cylinders of lengths a. This, for example, agrees with the Monte Carlo values for £(3i#3i) 
and e(3i#3t) to about 5%. 

Most of the knot energies in Table 1 have been taken from ^21, but we have independently 
calculated the energy of 2^, A\ and 6f exactly and the energy for several other knots and links 
approximately. We find £{2\) = An k, 12.57, to be compared with the Monte Carlo value 12.6. We 
also find e(4f ) = Q-k + 2 and £(6f ) = Stt -I- 3, where there are no Monte Carlo comparisons available, 
or needed. 

Model. — In our model, the chromoelectric fields Foi are confined to knotted tubes, each 
carrying one quantum of conserved flux |17j |18) . We consider a stationary Lagrangian density 

C ^ itrFo,i^« - y, (3) 

where, similar to the MIT bag model ^^Ij we included the possibility of a constant energy density 
V . To account for conservation of the flux $e, we add to J^jl the term 

trA{$£/(^a2) _ n'Fo,}, (4) 
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where is the normal vector to a section of the tube and A is a Lagrange multipUer. Varying the 
full Lagrangian with respect to , we find 

D''iFo,-Xn,)^0, D\Fo,-\n,)^Q, (5) 

which has the constant field 

Fo^ = ($ij/7ra2)n, (6) 

as its solution. With this solution, the energy is positive and proportional to I and thus the minimum 
of the energy is achieved by shortening Z, i.e. tightening the knot. 

We proceed to identify knotted and linked QCD flux tubes with glucballs, where we include all fj 
and f'j states. The lightest candidate is the /o(600), which we identify with the shortest knot/link, 
i.e., the 2\ link (sec Figure^; the /o(980) is identified with the next shortest knot, the 3i trefoil 
knot (see Figure and so forth. All knot and link energies have been calculated for states with 
energies less then 1680 McV. Above 1680 MeV the number of knots and links grows rapidly, and few 
of their energies have been calculated. However, we do find knot energies corresponding to known 
fj and f'j states, and so can make preliminary identifications in this region. (We focus on fj and 
f'j states from the PDG summary tables. The experimental errors are also quoted from the PDG. 
There are a number of additional states reported in the extended tables, but some of this data is 
either conflicting or inconclusive.) 

Our detailed results are collected in Tabled where we list fj and f'j masses, widths, and our 
identifications of these states with knots, together with the corresponding knot energies. 

In Figure|2|we compare the mass spectrum of / states with the identified knot and link energies. 
Since errors for the knot energies in were not reported, we conservatively assumed the error to 
be 1%. A least squares fit to the most reliable data (below 1680 MeV) gives 

i;(G) = (23.4 ±46.1) + (59.1 ±2. l)e(/-i:) [MeV], (7) 

with = 9.1. The data used in this fit is the first seven fj states (filled circles in Figure in the 
PDG summary tables. Inclusion of the remaining seven (non-excitation) states (unfilled circles in 
Figurc|2Jl in Tablc^ where either the glueball or knot energies are less reliable, does not significantly 
alter the fit and leads to 

i;(G) = (26.9 ± 24.9) + (58.9 ±1.0)e(if) [MeV], (8) 

with = 10.1. The fit (TJ is in good agreement [201 with our model, where E{G) is proportional to 
e{K). Better HEP data and the calculation of more knot energies will provide further tests of the 
model and improve the high mass identification. 

In terms of the bag model ^H], the interior of tight knots correspond to the interior of the 
bag. The flux through the knot is supported by current sheets on the bag boundary (surface of 
the tube). Knot complexity can be reduced (or increased) by unknotting (knotting) operations |21[ 
I22| . In terms of flux tubes, these moves are equivalent to reconnection events Hence, a 

metastable glueball may decay via reconnection. Once all topological charge is lost, metastability 
is lost, and the decay proceeds to completion. Two other glueball decay processes are: flux tube 
(string) breaking; this favors large decay widths for configurations with long fiux tube components; 
and quantum fluctuations that unlink flux tubes; this would tend to broaden states with short flux 
tube components. As yet wc arc not able to go beyond these qualitative observations, but hope to 
be able to do so in the future. 

We have assumed one fluxoid per tube. There may be states with more than one fluxoid, but these 
would presumably have somewhat fatter flux tubes with higher flux densities and higher energies. 
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Figure 1: The shortest knot/link soHtonic flux configuration has the topology of two linked tori, 
which in knot theory notation is 2f . This corresponds to the hghtest glueball candidate /o(600). 
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Figure 2: The second shortest soUtonic flux configuration is the trefoil knot 3i corresponding to the 
second hghtest glueball candidate /o(980). 
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Figure 3: Relationship between the glueball spectrum E{G) and knot energies s{K). Each point 
this figure represents a gluebaU identified with a knot or hnk. The straight fine is our model and 
drawn for the fit 0. 



For example, the two fluxoid trefoil knot 3i would certainly have e{K) > 2£:(3i) and a fairly reliable 
estimate gives s{K) « 2-\/2e(3i). Hence most multifluxoid states would be above the mass range of 
known glueballs. 

Discussions and conclusions. — In principle, lattice calculations can find any tame knot (knot 
without an infinite number of crossings or other pathology configuration, since there is always a 
contour through the lattice that represents the knotted path by some specific Wilson loop. However, 
since one is constrained by the rigidity of the lattice, energy minimization is difficult and requires 
a very fine-grained lattice. Thus we expect shape-evolving Monte Carlo techniques ^S] to be much 
more efficient and accurate for this purpose. 

Now we must discuss the details of identifications made in Tabled First, the /2(1270) does 
have a quark model interpretation. Either the glueball state in this range is well mixed with the 
quark model state and is part of the resonance at 1275 MeV (our interpretation), or the glueball 
state is yet to be discovered. In either case, more data in this region would be helpful. Next, the 
four (unconfirmed) glueball states with masses less then 1680 MeV from the extended PDG tables 
are identified as follows: (1) the 4f hnk with E{G) = 1289 and the 4i knot with E{G) = 1277 
are nearly degenerate, and the /i(1285) could actually be a pair of nearly degenerate states with 
identical quantum numbers associated with these knots; this is a possible interpretation of the 
/i(1285) mass measurements summarized on page 481 of Ref. 0; (2) the /2(1430) is treated as a 
rotational excitation of the /i(1420) and identified with the 5i knot; the energy difference between 
these two states, S', is a few MeV, but not well determined; this difference is of the order of what 
one would expect for rotational excitations; [We approximate E{fj) = E[fQ) -\- \J{J + 1)(5.] (3) 
we treat the /i(1510) as the ffist and the /2(1525) could be the second j^U rotational excitation 
of the /o(1500), which we identify with the 62 knot; now the energy step size is 5 5MeV which 
agrees with a simple estimate; (4) we assign the /2(1565) and the /2(1640) to the b\ and the 63 links 
respectively. 

Further details of knot excitations would be interesting to investigate, as would quantum and 
curvature corrections. At present we do not have a reliable way to estimate all these effects, nor 
do we have a good way to calculate glueball decays. However, we do expect high mass glueball 
production to be suppressed because more complicated non-trivial topological field configurations 
are statistically disfavored. 

Finally, knot solitons may also be able to survive within a quark-gluon plasma (e.g., in the interior 
of a RHIC event, quark star, or in the early universe). Complications will certainly arise in these 
cases due to additional parameters describing the media, as with knotted and linked electromagnetic 
plasma solitons; but if one holds the parameters constant throughout the region of interest in this or 
any sytem that supports knot/link solitons, the energy spectrum will be universal up to a scaling. 
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''Notation means a link of I components with n crossings, and occurring in the standard table of links (see e.g. 

1211 ') on the place. K#K' stands for the knot product (connected sum) of knots K and K' and K * K' is the link 
of the knots K and K' . 

''Values are from I15| except for our exact calculations of 2j, 2^ * Oi, and (2^ * Oi) * Oi in square brackets, our 
analytic estimates given in parentheses, and our rough estimates given in double parentheses. 
'^EiG) is obtained from e{K) using the fit 171 . 
•^States in braces are not in the PDG summary tables. 
'^This is the link product that is not 2^ * 2^. 

^Resonances have been seen in this region, but are unconfirmed [2]- 
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Tabic 1: Comparison between the glueball mass spectrum and knot energies. 



State 


Mass 


Width 


K ^ 


eiK) 


E{G) " 


/o(600) 


400 - 1200 


600 - 1000 


2? 


12.6 [in] 


768 [766] 


/o(980) 


980 ± 10 


40 - 100 


3i 


16.4 


993 


/2(1270) 


1275.4 ± 1.2 


185.1^^-^ 


2? *0i 


[6tt + 2] 


[1256] 


/i(1285) 


1281.9 ±0.6 


24.0 ± 1.2 


4i 
4? 


21.2 
(21.4) 


1277 
(1289) 


/i(1420) 


1426.3 ± 1.1 


55.5 ±2.9 


5i 


24.2 


1454 


{/2(1430) 


« 1430} 




5i 


24.2 


1454 + 6' 


/o(1370) 


1200- 1500 


200 - 500 


3i * Oi 


(24.7) 


(1484) 


/o(1500) 


1507 ±5 


109 ±7 


52 


24.9 


1496 


{./i(1510) 


1518 ±5 


73 ± 25} 


52 


24.9 


1496 + S 


/2(1525) 


1525 ± 5 


76± 10 


52 


24.9 


1496 + 36 


{/2(1565) 


1546 ± 12 


126 ± 12} 


5? 


(25.9) 


(1555) 


{/2(1640) 


1638 ±6 




6i 


((27.3)) 


((1638)) 



(22*0i)*0i<= [87r + 3] [1686] f 



/o(1710) 1713 ±6 125 ±10 6^ ((28.6)) ((1714)) 

3i#3t 28.9 (30.5) 1732 (1827) 

3i#3i 29.1 (30.5) 1744 (1827) 

2? * 2l [87r + 4] [1745] 

62 29.2 1750 
61 29.3 1756 

63 30.5 1827 
7i 30.9 1850 

819 31.0 1856 

820 32.7 1957 
/2(2010) "^QUtlo 202 ±60 72 33.2 1986 
/4(2050) 2025 ±8 194 ± 13 821 33.9 2028 

81 37.0 2211 

10i6i,i62 37.6 2247 

/2(2300) 2297 ±28 149 ± 40 818, 9i 38.3 2288 

/2(2340) 2339 ±60 3191?^ 92 40.0 2389 

lOi 44.8 2672 

111 47.0 2802 
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